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ABSTRACT 

Context. We present the results from a detailed X-ray variability analysis of 66 AGN in the Lockman Hole, which have 
optical spectroscopic identifications. 

Aims. We compare, quantitatively, their variability properties with the properties of local AGN, and we study the 
"variability - luminosity" relation as a function of redshift, and the "variability - redshift" relation in two luminosity 
bins. 

Methods. We use archival data from the last 10 XMM-Newton observations of the Lockman Hole field to extract light 
curves in the rest frame, 2-10 keV band. We use the "normalized excess variance" to quantify the variability amplitude. 
Using the latest results regarding the AGN power spectral shape and its dependence on black hole mass and accretion 
rate, we are able to compute model "variability - luminosity" curves, which we compare with the relations we observe. 
Results. When we consider all the sources in our sample, we find that their variability amplitude decreases with in- 
creasing redshift and luminosity. These global anti-correlations are less pronounced when we split the objects in various 
luminosity and redshift bins. We do not find a significant correlation between variability amplitude and spectral slope, 
r. The "variability - luminosity" relation that we detect has a larger amplitude when compared to that of local AGN. 
We also find that, at a given luminosity, the variability amplitude increases with redshift up to z~l, and then stays 
roughly constant. 

Conclusions. Our results imply that the AGN X-ray mechanism operates in the same way at all redshifts. The accretion 
rate (in units of the Eddington limit) for the objects in our sample increases from ~ 0.25, at z ~ 0.5, to 0.5 at z ~ 3. It 
does not exceed unity even in the case of the most luminous AGN. Their black hole mass also increases with redshift. 
The upper limit we find is consistent with the largest black hole masses found to date in the local Universe. The 
black hole mass increase, and the decrease of the rest frame light curve's duration with increasing redshift, can explain 
the global variability amplitude - redshift/luminosity anti-correlations that we observe. Among objects with the same 
luminosity, the black hole mass decreases and the accretion rate increases with larger redshift. This effect explains the 
increase of the variability amplitude up to redshift~l (at fixed luminosity bins). At higher redshifts, the decrease of 
the light curves length affects the variability amplitude significantly, forcing it to remain essentially constant. 
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1. Introduction 

The X-ray variability properties of nearby AGN have been 
extensively studied the last twenty years. One of the ear- 
liest, and best established so far, result was that the vari- 
ability amplitude decreases with increasing luminosity (e.g. 
Barr & Mushotzky 1986; Lawrence & Papadakis 1993; 
Green et al. 1993; Nandra et al. 1997; Turner et al. 1999; 
Markowitz & Edelson 2004) although at very low luminosi- 
ties the anti-correlation breaks, and the variability ampli- 
tude drops to nearly zero (e.g. Ptak et al. 1998). This "vari- 
ability - luminosity" anti-correlation may in fact be the re- 
sult of a more fundamental "variability - black hole (BH) 
mass" relation (e.g. Papadakis 2004; O'Neil et al. 2005). 
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If that is the case, then one could in principle get an esti- 
mate of the BH mass from a simple variability amplitude 
measurement (e.g. Nikolajuk, Papadakis & Czerny 2004; 
Nikolajuk et al. 2006). Both the "variability - luminosity" 
and the "variability - BH mass" anti-correlations are prob- 
ably projections of a "variability - BH mass - accretion 
rate" fundamental plane (M'^Hardy et al. 2006). 

The "variability BH mass estimation method" may be 
important in the case of distant AGN, which have been de- 
tected in recent X-ray surveys and in many cases they are 
too faint in other wavelengths for alternative BH mass es- 
timation methods to apply. However, until recently, knowl- 
edge of the X-ray variability properties of distant (z>0.1) 
AGN was scarce. The situation changed the last few years, 
when data from deep X-ray surveys become available, and 
techniques that are able to measure variability amplitudes 
for low signal-to-noise sources were developed. 
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Almaini et al. (2000) and Maimers ct al. (2002) pre- 
sented the results from an analysis of the X-ray variabil- 
ity of radio-quiet AGN selected from a deep ROSAT sur- 
vey and the i?05'ylT' archive, respectively. They confirmed 
the decline in variability amplitude with luminosity out to 
redshifts as large as z~2. They also noted some evidence 
for an increase in the variability amplitude at high red- 
shifts. Paolillo et al. (2004) studied the X-ray variability of 
sources detected in the Chandra Deep Field South. They 
also observed an anti-correlation with the luminosity of the 
sources and, in agreement with the previous studies, they 
found that high-redshift objects have larger variability am- 
plitudes than expected from extrapolation of their low-z 
counterparts. 

In the current study we measure the variability ampli- 
tude of AGN detected by XMM-Newton in the Lockman 
Hole (LH) field. XMM-Newton has carried out its deepest 
survey so far in the direction of this field. Mateos et al. 
(2005) studied the spectral properties of the 123 bright- 
est sources (i.e. sources more than 500 MOS-I-PN counts 
in the 0.2-12 keV band) in this field. Recently, Mateos et 
al. (2007) presented the results from a detailed study of 
their flux variability properties as well, using data from 16 
of the currently available XMM-Newton observations. They 
found that the variability amplitude does not significantly 
depend on redshift or X-ray luminosity. They also studied 
their spectral variability properties, and found that spec- 
tral variability is less common that flux variability and a 
lack of correlation between flux and spectral variability. 

In this work we study the variability properties of 66 
objects with optical spectroscopic identifications in the LH 
field, using data from the last 10 XMM-Newton observa- 
tions of the LH field only. These observations were per- 
formed over a period of two months on a quasi-regular pat- 
tern. As a result, we were able to construct uniformly sam- 
pled, high signal-to-noise light curves for all sources in our 
sample. Furthermore, different to all previous studies, we 
extracted rest frame 2-10 keV band light curves, and we 
compare, quantitatively, our results with those from the lo- 
cal AGN sample of Markowitz & Edelson (2004). We also 
quantify the variability amplitude in a way different to what 
has been used in most cases in the past (i.e. Almaini et al. 
2000, Manners et al. 2002, Mateos et al. 2007): we use the 
"normalized excess variance" a rather well known and fre- 
quently used estimator which is easy to interpret. 

Our main aim is to study the correlation between 
the variability amplitude and other source parame- 
ters like redshift, X ray luminosity and X ray spectral 
slope and, in particular, to investigate if and how the 
"variability-luminosity" and "variability - redshift" corre- 
lations changes with redshift and luminosity, respectively. 
Our results suggest that the X-ray source in AGN operates 
in the same way at all redshifts, and allow us to estimate 
how the average black hole mass and accretion rate of the 
AGN in our sample change with redshift. 



2. The XMM-Newton observations and data 
reduction 

XMM-Newton observed the Lockman Hole field 17 times 
from April 2000 to December 2002. The total number of X- 

ray sources detected in this field is 268 (Mateos et al., 2005). 
Currently, 74 of these objects have been spectroscopically 



Table 1. Summary of the 10 LH XMM-Newton observa- 
tions that we used in our study. In the last column we list 
the "good time interval" (GTI) after the removal of high 
background periods. 



ObsID 


Obs. 


Datc 


GTI (ksec) 


0147510101 


2002- 


-10- 


-15 


91 


0147510801 


2002- 


-10- 


-17 


75 


0147510901 


2002- 


-10- 


-19 


90 


0147511001 


2002- 


-10- 


-21 


82 


0147511101 


2002- 


-10- 


-23 


50 


0147511201 


2002- 


-10- 


-25 


104 


0147511301 


2002- 


-10- 


-27 


65 


0147511601 


2002- 


-11- 


-27 


100 


0147511701 


2002- 


-12- 


-04 


99 


0147511801 


2002- 


-12- 


-06 


88 



identified as either Type I or Type II AGN. They constitute 
our initial sample of sources. In order to investigate their 
variability properties we used data from the last 10 XMM- 
Newton observations only. A summary of these observations 
is reported in Table 1. 

By choosing to work with observations which are spread 
over a period less than two months we minimize the risk 
that potential changes in the response of the XMM-Newton 
detectors over large period of time will add to the observed 
variations in the light curves. Furthermore these 10 ob- 
servations were done in a quasi-regular temporal pattern. 
The first seven were performed every 2 days over a period 
of 14 days while the last three over a period of 10 days. 
Consequently, the choice to use data from these observa- 
tions only, minimizes differences in the sampling pattern of 
the resulting light curves. 

2.1. Data reduction 

During the XMM-Newton observations the EPIC PN and 
MOB cameras were operated in full frame mode with 
a medium filter. The EPIC data were reprocessed with 
XMMSAS version 6.5. We have selected PN photons with 
PATTERN< 4 (i.e. singles and doubles) and FLAG=0. In 
the case of the MOS cameras we have selected events with 
PATTERN < 12 and FLAG=0. 

Source counts were accumulated from a circular region 
of 18" radius, centered around the coordinates of each 
source (as listed in Mateos et al. 2005), in the energy range 
from 2 to 10 keV, in each source's rest frame. The back- 
ground count rate, for all EPIC detectors, was determined 
by accumulating counts from six circular, source free re- 
gions, each one of radius 85". Data from periods of high 
background were disregarded from further analysis. Given 
the different redshifts, we constructed the appropriate 2 — 10 
keV background light curves for each source individually. 
Finally, we produced PN, MOSl and M0S2, background 
subtracted, rest frame 2 — 10 keV light curves, for each 
source in our sample. 

Since the XMM-Newton EPIC cameras have different 
geometry, it was common that a number of sources detected 
in each observation would fall near or inside CCD gaps, 
and/or close to the edge of at least one of the detectors. 
Furthermore, in some cases we would detect bad/hot pix- 
els within a source's photon extraction aperture in one (or 
more) of the EPIC cameras. For these reasons, we inspected 
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Table 2. IDs, redshift, spectral and variability properties of the 66 AGN in the Lockman Hole field that we study in 
this work. The first column lists the source ID as given by Mateos et al. (2005) together with its RA and Dec. Data 
listed in columns 2, 9, and 10 are taken from Table 8 of Mateos et al. (2005). Numbers in parenthesis in the third column 
correspond to the 3cr upper limits on cr^xS' ^^^^ '-'^ NV-sources. The values listed in column 9 correspond 

to the logarithm of the rest frame, 2-10 keV band luminosity of each object. Column 8 indicates whether a light curve 
is considered as "variable" (V) or "non- variable" (NV) (for details see Section 3.1). Column 11 lists the XMM-Newton 
detectors that were used to construct the respective light curve. The full table is available in electronic form at the CDS. 



SourcelD 


z 


'2 

"NXS 


log o-Nxs 


xVdof 


V/NV 


log(Lx) 


r 


Instr. Used 


5(10 52 43.3 +57 15 45.9) 


2.144 


0.026+0.015 


-1.59+0.26 


7.2/4 


NV 


44.68 


1.90 


PN 




(0.045) 












39(10 53 19.09 +57 18 53.6) 


0.711 


-0.019+0.027 




1.8/4 


NV 


43.44 


1.79 


M1+M2 






(0.081) 














41(10 51 19.14 +57 18 34.1) 


1.640 


0.445+0.089 


-0.35+0.09 


25.5/5 


V 


44.11 


2.06 


PN 



visually all sources in each PN, MOSl and M0S2 image. 
We would also check whether the background count rate we 
had estimated in an automatic way would be representative 
of the local background around each source. 

As a result of this screening process, light curves from 
each EPIC detector have a different number of points in 
most cases. Due to the geometry of the PN detector, the 
number of sources falling on the PN CCD gaps is larger 
than the MOS cameras. Consequently, the number of data 
points, Ndata, in the PN light curves is usually smaller than 
the number of points in the MOS light curves. 

A reliable estimate of the variability amplitude requires 
the use of light curves with: a) the largest possible Ndata 
and b) the highest signal-to-noise ratio. This ratio can be in- 
creased by combining counts from all three EPIC camera^. 
We therefore added the background subtracted PN and 
MOS light curves to create combined PN+M0S1+M0S2, 
PN+MOSl, PN+M0S2 and M0S1+M0S2 light curves. 
Among the individual and combined light curves of each 
source we would finally choose the one with the largest 
number of points. In the case there were more than one, we 
would choose the light curve with data taken from as many 
as possible EPIC detectors. 

The largest Ndata value among all light curves is eight. 
Since Ndata (and hence the sampling pattern as well) should 
be roughly the same in all of them, we required that 
Ndata > 4 (i.e. half of Ndata.max)- There are 66 (out of the 
74) AGN whose light curves satisfy this criterion. These 
sources constitute our final sample. Table 2 lists their ID 
number, redshift. X-ray spectrum slope, P, and rest frame, 
2-10 keV band luminosity, Lx. In the last column of the 
same Table we list the combination of the EPIC cameras 
that we used for each source. 

There is just one object with Ndata = 4 in our sample, 
and four with Ndata — 8. In most cases, Ndata =6. Due to 
this fact, the sampling pattern of most light curves is in- 
deed similar, in the observer's frame. However, when one 



^ One must be careful in combining data from the three EPIC 
cameras, given the difference in their response. However, the 
uncertainty introduced to our results by the difference in the 
average energy of the (rest frame) 2-10 keV band photons in 
the EPIC PN, MOSl and M0S2 cameras is smaller than that 
introduced by say the small number of points in the light curves 
and the (unknown) stochastic nature of the underlying variabil- 
ity process. For this reason, we use light curves with photons 
accumulated from as many cameras as possible, as explained in 
the text. 



takes into account the different redshifts of the sources in 
our sample, the light curves' duration changes, from one 
source to the other, in their rest frame. This redshift in- 
duced inhomogeneity in the temporal pattern of the light 
curves has to be taken into account when we compare and 
interpret the variability properties of the low and high-z 
objects in our sample. 

3. Variability analysis method 

3.1. Detection of intrinsic flux variations 

First we investigated which light curves show significant 
variations, using the test. Since our sample consists 
of the brightest sources in the catalogue of Mateos et al. 
(2005), there are more than 15 source photons in each point 
of their light curves. This guarantees the applicability of 
Gaussian statistics, and hence the reliability of the test 
results. The results (x^ /degrees of freedom) are listed in the 
7th column of Table 2. We accept that a source shows sig- 
nificant fiux variations if the probability of being larger 
than the obtained value, by chance, is smaller than 5%. 
Column 8 on the same table lists whether a source is vari- 
able ( "V" ) or not ( "NV" ) , according to this criterion 0. 

Sixteen objects, i.e. 24% of the sources in the sample, 
turned out to be non variable. In order to investigate this 
issue further, we used the Kolmogorov-Smirnov (K~S) test 
to compare the distribution functions of the redshift, Ndata 
and mean count rate (CR) of the V and NV-subsamples. 

We find that the redshift distribution functions of the 
two subsamples are almost identical. However, we get a 
different view when we consider the results from the com- 
parison of the CR and Ndata distributions. The NV-sources 
are systematically fainter than the V-sources. We find that 

< CRnv >= 7.1 X 10"'* cnts/s and < CRy >= 2.2 x lO^^ 
cnts/s. The K-S test suggests that the difference in bright- 
ness is significant at the 96% level. There are also less 
points in the light curves of the NV-sources. We find that 

< Ndata,NV >= 5.7 as opposed to < Ndata,v >= 6.2. The 
K-S test suggests that the difference in the distribution of 
Ndata is significant at the 99.94% level. We conclude that 
we do not detect intrinsic variations in the NV-sources most 
probably because they are less bright and their light curves 
have fewer points than the V-sources. 

^ In the Appendix we discuss in detail the possibility that the 
observed variations may be affected significantly by instrumental 
effects 



^ 1.5 

U 

1 1 

c 

S 0.5 

X 

2 1-5 
? 1 
5 0.5 



V364 



V321 



10 15 
Time (days) 



25 



Fig. 1. The rest frame, 2-10 keV light curve of source V364, 
which is a "NV" object, and V321, which belongs in the "V" 
subsample. Time is measured in days (in rest frame) since 
the first observation of each source. 

Fig. 1 shows the rest frame, 2-10 keV light curve of the 
objects with ID numbers V364 and V321. Their average 
count rate is similar, and in fact there are more points in 
the V364 light curve. Nevertheless, it is source V321 which 
shows significant variations in its light curve. As Fig. 1 
shows this is due to just one point (which we have carefully 
examined to establish its reliability). This result demon- 
strates clearly the importance of working with light curves 
which have the largest possible number of points (and jus- 
tifies the way we constructed the light curves, as explained 
in the previous section). 



3.2. Determination of variability amplitude 



We use the normalized excess variance, u^xs (^-S- Nandra 
et al. 1997) as a measure of the intrinsic variability ampli- 
tude of the light curves. This estimator is defined as. 



'NXS 



(1) 



where < x >, and 5^ = (l/Njata) Ei=T"''(xi- < x >)^ are 
the mean and variance of the light curve, while < cr^j.^. >= 

(1/Ndata) Er=r ^err i average contribution of the 

Poisson noise process to the observed scatter around the 
mean. Its square root (i^var) indicates the average variabil- 
ity amplitude of a source as a fraction of its light curve 
mean. 

The spread of the observed cr^xs values around the in- 
trinsic, "true" tJ^ value (i.e. the "error" of cr^xs) should 
decrease with increasing Njata- This spread also depends 
on the stochastic nature of the process underlying the AGN 
X-ray variability (see Vaughan et al. 2003 for a detailed dis- 
cussion on this issue). In fact, even if a source is not intrin- 
sically variable, the cr^xs estimates based on light curves 
from observations performed at different times will not be 
identical, due to the presence of the Poisson noise alone. 

It is not easy to estimate a priori how does the error on 
"'nxs depend on Njata and on the stochastic nature of the 
variability mechanism. On the other hand, the error due to 
Poisson noise can be estimated using the following equation 
(Vaughan et al. 2003), 



eiT(cr^xs) 



N 



data 



N 



data 



< X 



(2) 
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Fig. 2. Plot of ctnxs ^ ^ function of redshift (top panel), 
logarithm of the rest-frame 2-10 keV luminosity (middle 
panel) and spectral slope value (bottom panel), for the 50 
variable sources in our sample. The vertical gray lines in 
the top panel show the 3(7 upper limits on tr^xs case 
of the NV-sources. 



This is useful in order to estimate upper limits on the in- 
trinsic excess variance in the case when we do not detect 
significant variations. 

Columns 3 and 4 in Table 2 list the excess variance, 
and its error, for all the sources in our sample. In the 
case of the NV-sources, Unxs consistent with zero, as 
expected. For these sources, the numbers in parenthesis in 
the third column correspond to the 3a upper limits, i.e. 
^Nxs + 3err((T^xs)- columns 5 and 6 we hst the log- 
arithm of cr^xs ^"^d its error, estimated using the usual 
"propagation of errors" formula. 



4. Correlation of variability amplitude with source 
parameters 

We now proceed to study the correlation between the ex- 
cess variance and source parameters such as the redshift, 
Lx, and F. In order to quantify the correlation between two 
parameters, we use the Kendall's t rank correlation coeffi- 
cient. We accept a correlation to be significant if the null 
hypothesis probability, Pnuii, is less than 5%. 

Filled squares in the upper panel of Fig. 2 show the 
""'nxs vs z" plot for the 50 variable sources in our sample. 
Errors are estimated using equation (2). A trend of de- 
creasing variability amplitude with increasing redshift can 
be seen. When we consider the [log(crNxs)' log(2)] data for 
the variable sources only, we find that r — —0.2, a result 
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which confirms the significance of the "variabihty - z" anti- 
correlation (PnuU = 0.04). 

The vertical grey lines in the same panel of Fig. 2 show 
the 3a upper limits on cr^xs ^^"^ '^'^^'-^ '-'^ NV-sourccs. 
They arc so large that they do not impose any significant 
constrains on the "cr^xs ~ ^" relation. We need better light 
curves for these objects (in terms of signal-to-noise ratio 
and Ndata) to constrain this relation better. The same holds 
true for all relations we study below. For that reason, we 
do not consider hereafter the NV-sources in our analysis. 

The solid line in the same figure shows the best fitting, 
"ordinary least squares bisector" line to the data in log-log 
space (the best fitting parameter values were estimated as 
in Isobe ct al. 1990). We find that a relation of the form 
"'nxs °^ ~ —1.05 ± 0.10 fits the data well. 

Next, we examined the relation between variability am- 
plitude and Lx. The middle panel in Fig. 2 shows the "o"xxs 
vs log(Lx)" plot for the V-objects in our sample. Given the 
correlation between Lx and redshift, it is not surprising 
that we an anti-correlation between variability amplitude 
and luminosity can be seen. Kendall's r suggests that the 
anti-correlation between log(CT^xs) ^^'^ log(Lx) is signifi- 
cant (r = -0.22, P„uu = 0.023). The sohd line shows the 
best fitting a^xs *^ ^x " line (a = —0.66 ± 0.12), which 
fits the data rather well. 

In the bottom plot we show the "cr^xs vs F" rela- 
tion for the V-sources in our sample. One could perhaps 
suggest a positive correlation between the two quantities, 
but Kendaull's r in this case does not support this claim 
(r = 0.03, P„uu = 0.74). The best fitting cr^xs « 
model to the data (a = 1.58 ± 1.14) is also shown with the 
solid line. The uncertainty is so large that the possibility of 
a = (i.e. of no correlation between cr^xs ^) cannot 
be excluded. 

The "ct^xs " log(Lx)" anti-correlation that we find is 
in agreement with what has been observed in nearby AGN, 
both on long time scales (e.g. Papadakis 2004; Markowitz 
& Edclson, 2004) and short time scales (e.g. O' Neil et 
al. 2005, and references therein). In fact our results indi- 
cate that the decline in variability amplitude with lumi- 
nosity holds out to redshifts as high as ~ 3. This decline 
can also explain the "t^xs ~ ^" anti-correlation: as the 
more distant objects are systematically more luminous, it 
is not surprising that they are also less variable. However, 
it is important to note here that, although the correlations 
we detect are statistically significant, they are also driven 
by a few individual points. For example, if we omit the 
high-z source with cr^xs ^ '-'•01 then Pnuii > 0.05 in the 
case of the [log(crNxs)5 log(z)] anti-correlation. We reach 
the same conclusion if we omit the two pints with the 
smallest ct^xs i^i the case of the [log(cr^xs)' log(Lx)] anti- 
correlation. Consequently, one has to be careful regarding 
the reality of these correlations. 

5. Comparison with nearby AGN 

In this section we present the results from a quantitative 
comparison between the "variability - luminosity" rela- 
tion as determined for nearby AGN and the V-objects in 
our sample, grouped in various redshift bins. We use the 
Student's-t test to compare the mean of two sample dis- 
tributions. We consider two means as being significantly 
different if PnuU is smaller than 5%. 



In order to determine the "log(<J^xs) log(Lx)" rela- 
tion for the AGN in the local universe we considered the 
12 objects in the sample of Markowitz & Edelson (2004). 
Object names and the dates of the .RXTi? observations that 
were used to extract light curves are listed in Table 1 of 
Markowitz & Edelson (2004). Details of the data reduction 
can also be found in the same paper. 

Filled squares in all panels of Fig. 3 show the [log(cr^xs)' 
log(Lx)] data of these AGN. Excess variances were esti- 
mated as described in Section 3.2. In the case when there 
were more than one light curves for an object, we plot 
the straight mc^an of the individual cr^xs values. We used 
1.6 days binned, 2 12 keV light curves, which were 23.2 
days long. This is equal to the average, rest frame light 
curve length of the V-sources in our sample. As for Lx, 
we adopted the 2 — 12 keV luminosity measurements of 
Markowitz & Edelson (2004). This is slightly different to 
the energy band we chose in this work, but this should not 
affect our results significantly. 

The open circles in the top, middle and bottom panel 
of Fig. 3 show the [log(CT^xs)! iog(Lx)] data of the V- 
sourccs with 0.2<z<l (the "Low" -redshift bin; 18 objects), 
l<z<2 ("High" -redshift bin; 21 objects) and 2.3<z<3.4 
("very High", or vHigh, redshift bin; 9 objects). The av- 
erage redshift of the objects in the Low-z, High-z and v- 
High bins is < zlow-z >= 0.7, < znigh-z >= 1-48 and 

< ZvHigh-z >= 2.85. 

In order to reduce the scatter of the Low-z data points 
in Fig. 3, we sorted them in order of increasing luminos- 
ity, we considered three bins of 6 points each, and we esti- 
mated their average log(crNxs) ^rid log(Lx) values (shown 
with filled circles in the top panel of Fig. 3). Similarly, filled 
circles in the middle panel show the mean log(cr^xs) and 
log(Lx) values of the data in three bins (each of 7 points) 
into which we grouped the High-z objects (in order of in- 
creasing luminosity). We also binned the vHigh-z objects in 
two groups and we estimated their average log((J^xs) ^^'^ 
log(Lx) values. The results are shown with the filled circles 
in the bottom panel of Fig. 3. 

On average, Low-z and nearby AGN have similar X- ray 
luminosity. We find that < log(Lx,Low-z) >= 43.62 as op- 
posed to < l0g(Lx,„earby) >= 43.31 (PnuU = 0.21). On the 
other hand, the High-z and vHigh-z AGN are significantly 
more luminous than the nearby AGN: < log(Lx, High-z) >= 
44.24 and < log(Lx,vHigh-z) 44.63 (P„uii = 4 x lO'^, 
and PnuU = 7.3 x 10~^, respectively). If both the LH and 
the nearby AGN followed the same "c^xs ^ log(Lx)" rela- 
tion, we would expect the Low-z objects to be as variable 
as the nearby AGN, and the High-z and vHigh-z AGN to 
be significantly less variable. However, this is not the case. 

The average variability amplitude of the nearby and 
the Low-z AGN are < log(cr^xs,nearby) —1.45 and 

< log(o-^xs,Low-z) >= -0.97. We find that their dif- 
ference is significant (Pnuii = 0.003). Even the High-z 
sources are significantly more variable than the nearby 
AGN (< log(a2xs,High-z) >= -1-13, PnuU = 0.031), while 
the vHigh-z objects are at least as variable as the nearby 
AGN, despite the fact that they are much more luminous 
(< log(a2xs,,High-z) >= -1-34, PnuU = 0.57). 

We conclude that, at a given X-ray luminosity, the 
z>0.2, V-sources in our sample are systematically more 
variable than the nearby AGN. This is in agreement with 
the results of Almaini et al. (2002), who found that the 
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Fig. 3. The "variability - luminosity" plot for the nearby 
AGN (filled squares in all plots), the Low-z, High-z and 
vHigh-z objects in the V-subsample (open circles in the 
top, middle and bottoms panels, respectively). Filled, grey 
circles correspond to the mean variability amplitude and 
luminosity, when data are grouped in various bins as ex- 
plained in the text. The solid line in the top panel and the 
long-dashed lines in all panels show the model "variabil- 
ity amplitude - log(Lx)" lines, estimated as explained in 
Section 6. 

z> 0.5 AGN do not show the anti-correlation with lumi- 
nosity seen in local AGN, and the results of Paolillo et al. 
(2004), who found that high-z objects have larger variabil- 
ity amplitudes than expected from extrapolations of their 
low-z counterparts. 

6. Impact of BH mass and accretion rate to the 
variability - luminosity relation 

Let us assume that the X-ray variability mechanism is the 
same in all AGN, irrespective of their redshift. This as- 
sumption implies that the power spectral density, PSD, has 
the same shape in all of them. Recent work, based on de- 
tailed analysis of high quality RXTE and XMM-Newton 
light curves, has shown that the AGN PSDs have a power- 
law shape of slope ~ —2 for the 2-10 keV band (as op- 
posed to ^ —2.7 for softer bands), which changes to a 
slope of ~ — 1 below a so-called "break frequency", I'tj 
(e.g. Markowitz et al., 2003; M^Hardy et al., 2004). 

The excess variance estimated from a light curve of 
length T, is an approximate estimator of the following in- 
tegral, 

POO 

4xs = / PSD{,y)d,^, (3) 

Ju,f 



where i^if — 1/T. If, as we assumed above, PSD{v) — 
A{v/vbf)~'^ {v < Vbf) and PSD{v) = A{v/vbf)^'^ (when 
V > Vbf), then, 

o'nxs = PSDa,np^^{i'bf) ~ ln(j^;/) + 1], (4) 

where PS Damp = Ai^bf 0.02 (Papadakis 2004). In the 
case when vif > i/bf, 

C^NXS = {PSD)arnp{vbflvif)- (5) 

Recently, M'^Hardy et al. (2006) have demonstrated that 
Ubf depends on both the BH mass, Mbh, and the accretion 
rate, riiEdd (measured in Eddington units), as follows, 

Ubf = 0.02977mEdd(MBH/10^MQ)-i, (6) 

where "q is the efficiency of the mass to energy conversion 
(hereafter we assume 77 — 0.1). One can now use the above 
equation to substitute Vbf in equations 4 and 5 in order 
to relate cr^xs "^ith Mbh and riiEdd- The resulting relation 
implies that for a given Mbh, ctnxs increases with increas- 
ing accretion rate, while, among all AGN with the same 
riiEdd, objects with smaller Mbh should have larger vari- 
ability amplitude. 

The cr^xs ~ /(^^BH, riiEdd) relation can be trans- 
formed to a cr^xs ~ /(Lx) relation if we use the fact 
that the bolometric luminosity, Lboi, is given by Lboi = 
1.3rymEddlO'^^(Mbh/Mo) erg/s, and we adopt an X-ray lu- 
minosity to Lboi conversion factor. To this end, we adopted 
the Marconi et al. (2004) prescription, 

log(Lx) = log(L*) - 1.54 - Q.2AI - 0.012^2 _ o.0015£^ (7) 

where £ = logL* — 12, L* being the bolometric luminosity 
in solar luminosity units. 

To summarize, it is possible to derive a model "c^xs 
Lx" relation, if we assume that: 1) all AGN, at all redshifts, 
vary like the local AGN, 2) the efficiency is the same in all 
objects, and 3) a bolometric to X-ray luminosity conver- 
sion relation. In addition, the model takes into account the 
different rest frame length of the derived light curves ( "the 
redshift induced inhomogeneity" we mentioned in Section 
2). Although Tobs is similar in all light curves, vif in equa- 
tions (4) and (5) should be equal to l/Ti-cstframc, where 

^rcstframc — Pohs/ (i- ^ ^)- 

The solid line in the upper panel of Fig. 3 shows such 
a model "log(cr^xs) log(Lx)" curve in the case when 
riiEdd = 0.05 and i^i/,nearby = 1/23.2 days~^. Clearly, this 
line describes rather well the "variability - luminosity" re- 
lation for the nearby AGN. This is not surprising given the 
fact that the model PSD shape that we have assumed, and 
the "i^bf vs BH mass and accretion rate" relation that we 
use, have resulted from the study of objects which are all 
included in our nearby AGN sample. 

The long dashed line in the same panel show the model 
"log(tT^-j^g) vs log(Lx)" relation in the case when riiEdd = 
0.25 and i^if — 1/29.6 days^^ (29.6 days is the average, rest 
frame duration of the Low-z light curves). This is entirely 
consistent with the mean [log(fT^xs)' ^°s(Lx)] data of the 
Low-z objects. 

The model curves plotted in the upper panel of Fig. 3 
can be used to explain the differences we observe between 
the nearby and Low-z AGN. For example, let us consider 
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Fig. 4. The "variability - redshift" relation for the variable 
sources in our sample with "low" and "high" luminosity 
(open circles in the top and bottom panels, respectively). 
Filled squares show the model predictions, estimated as ex- 
plained in Section 6.1, and stars indicate the average vari- 
ability amplitude and redshift of the "lowLum" and "high- 
Lum" , nearby AGN. 



the objects with log(Lx) — 44. According to the model, in 
the local Universe these are objects with log(MBH) ~ 8.7 
and rii ~ 5% of the Eddington limit. The fact that the re- 
spective Low-z sources are more variable can be explained 
by a smaller black hole mass (log(MBH) ~ 8) and a higher 
accretion rate (rh ~ 25% of the Eddington limit). At the 
same time, the higher accretion rate compensates the dif- 
ference in BH mass, and explains why objects with different 
BH mass have the same luminosity. 

The long dashed lines in the middle and bottom panels 
of Fig. 3 show the model "log(crNxs) '^^ log(Lx)" lines for 

rilEdd.High-z = 0.4, hlEdd.vHigh-z = 0.5, ly'lf Migh-z = 1/20.2 

days"-^ and i'//,vHigh-z — 1/12.7 days^^ (20.2 and 12.7 
days being the average, rest-frame length of the High-z 
and vHigh-z light curves, respectively). Clearly, these model 
curves describe well the mean variability amplitude vs X- 
ray luminosity data of these objects. 



6.1. The "variability - redsliift" relation in AGN 

According to the model presented above, as the redshift in- 
creases, objects at a given luminosity bin should correspond 
to AGN of smaller BH mass (hence the larger variability 
amplitude) and higher accretion rate. Therefore, we would 
expect the variability amplitude of AGN with the same 
luminosity to increase with increasing redshift. In order to 
investigate this issue further, we considered the "variability 
- redshift" relation at fixed luminosity bins. 

Open circles in Fig. 4 shows the [log((T^-xs)'l°g(^)] 
data for the variable objects in our sample with Lx = 
X0''3-6 _ 10441 ergs/s (top panel; the "lowLum" bin) and 
Lx = lO'^'*-! - lO^''-^ ergs/s (bottom panel; the "highLum" 



bin). Star-like points in both panels show the average vari- 
ability amplitude and redshift for the nearby, "lowLum" 
and "highLum" AGN. Contrary to the 'variability - red- 
shift" anti-correlation we detect when we consider all vari- 
able sources together (see top panel in Fig. 2), here one can 
even claim a positive correlation between cr^xs ^^'^ 
shift: the variability amplitude of the higher redshift, LH 
AGN is higher than the amplitude of the nearby objects 
with the same luminosity. This is in agreement with the 
discussion in the paragraph above. 

However, the situation is less clear when we consider 
the LH sources only (i.e. the AGN at redshift larger than 
^ 1 or so). We find that there does not exist a signifi- 
cant anti-correlation (or correlation) either in the "low" or 



"highLum" bins (riowL u 



-0.16, ThighLu 



-0.17, with 



Pnuii.iowLum = 0.4, and Pnuii.highLum = 0.37, respectively). 

The average X~ray luminosity of the lowLum and 
highLum sources is < log(Lx,iowLum) >= 43.8 and < 
log(Lx, highLum) >= 44.3. We can now use the riiEdd — 
0.25, 0.4 and 0.5 model curves, shown with the dashed lines 
in Fig. 3, to estimate the amplitude of objects with Lx — 
IQii.s g^j^j _ lO'*'' '^ ergs/s in the Low-z, High-z and 
vHigh-z redshift bins, respectively. Filled (grey) squares 
in Fig. 4 show the variability amplitude of these objects 
(as X-axis coordinates we have used the average redshift of 
all sources in the Low-z, High-z and vHigh-z bins, while 
the horizontal bars correspond to the size of the respective 
bins). Clearly, the model "variability - redshift" relations 
agree well with the data in both luminosity bins. 

The main reason that we do not observe a strong, posi- 
tive correlation between variability amplitude and redshift 
(when z>l) at fixed luminosity bins is that, the rest frame 
light curve length, Trostframc, and hence o-^xs well, de- 
creases with increasing redshift. Therefore, although the 
more distant AGN have a smaller BH mass and higher ac- 
cretion rate than the local, same luminosity objects, their 
variability amplitude is not as large as we would expect 
due to this decrease in Trcstframo- This "Trostframc - red- 
shift" anti-correlation is strong enough to compensate the 
increase in cr^xs '^^'^ ^'^ ^^^'^ higher riiEdd and smaller Mbh- 



7. Discussion and conclusions 

We have studied the long term X-ray variability properties 
of 66 AGN in the Lockman Hole field which have optical 
spectroscopic identifications. We used rest frame 2-10 keV 
band light curves, and we have taken particular care that 
they all have a similar sampling pattern. For this reason, we 
extracted data from the last 10 XMM-Newton observations 
of this region only, and we restricted our sample to objects 

with Ndata > 4. 

Using the test, we detected significant variations in 
50 of them. The low signal-to-noise ratio and the smaller 
number of points in the light curves of the remaining 16 
objects can explain the absence of significant variations in 
them. We have used the well known "normalized excess 
variance" , (j'^xs ' estimator in order to quantify the vari- 
ability amplitude of these AGN. We have also used RXTE 
light curves of 12 nearby AGN, with an average energy and 
length similar to the LH AGN, in order to measure their 
normalized excess variance and compare their variability 
properties with the Lockman Hole AGN. Our results can 
be summarized as follows: 
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a) When wo consider all the objects in our sample 
together, we detect a significant anti-correlation between 
redshift, luminosity and variability amplitude: the ampli- 
tude decreases with increasing redshift and luminosity. This 
"variability luminosity" anti-correlation is in agreement 
with (and extends out to redshift of ~ 3) the results of 
Almaini et al. (2000), Manners et al. (2002) and PaoUUo et 
al. (2004). 

b) We do not find a significant correlation between vari- 
ability amplitude and F. This is in agreement with the re- 
sults of Mateos et al. (2007) but opposite to what has been 
observed in nearby AGN, where sources with steeper spec- 
tra show larger amplitude variations (e.g. Green et al. 1993, 
Turner et al. 1999, Grupe ot al. 2001). Paolillo ct al. (2004) 
also observed sources with harder energy spectra to be less 
variable in their study of the X-ray variability properties 
of the AGN in the Chandra Deep Field-South. In order to 
investigate this issue further, we searched for a "variability 
- F" relation in the Low-z, High-z and vHigh-z subsamples, 
but we did not detect any. Perhaps, if real, such a "variabil- 
ity - spectral shape" correlation has such a large intrinsic 
scatter, in which case we would need a larger sample to 
detect it. 

c) The "variability - luminosity" relation of the 
Lockman Hole AGN has a larger amplitude when compared 
to that of local AGN. This is in agreement with the results 
of Almaini ct al. (2000), Manners ct al. (2002) and Paolillo 
et al. (2004), who also found that high-redshift sources 
behave differently, in terms of their variability amplitude, 
than their nearby counterparts. 

d) The global "variability amplitude- 
redshift/luminosity" anti-correlations arc loss pronomiccd 
when we consider AGN in various luminosity and redshift 
bins. The main trend we detect is that, at a given luminos- 
ity bin, the variability amplitude increases with redshift 
until Z'^l, and then stays roughly constant. 

It is not possible to compare rjuantitatively our results 
with those presented in the past. Quantitative comparisons 
can be performed when one uses light curves of similar dura- 
tion and rest frame energy band. Furthermore, the variabil- 
ity amplitude must be measured in the same way, prefer- 
ably with the use of an estimator whose relation with the 
AGN fundamental physical parameters, like the BH mass 
and accretion rate, can be established. Such an estimator 
is the "normalized excess variance" , which wc adopted in 
this work, since its dependence on Mbh, hiEdd (^nd the light 
curve's length can be easily determined (Section 6.1). 

In any case, we believe that parameters like the slope 
of the best fitting line in the overall "variability- lumi- 
nosity" or "variability redshift" plots (top and middle 
panels in Fig. 2) do not provide important physical in- 
sights. The shape of the "variability-luminosity/redshift" 
relations is strongly affected by factors like the length of 
the light curves used and their energy band. Such "exper- 
imental" parameters can introduce or even "destroy" any 
global "variability- luminosity" or "variability - redshift" 
anti-correlations. For example, as we discussed in Section 4, 
although formally significant, the "variability luminosity" 
and "variability - redshift" anti-correlations that wc detect 
are mainly driven by a few points. A better understanding 
of how the AGN variability properties change with look- 
back time can be achieved by studying the "variability- 
luminosity" and/or the "variability redshift" relations in 
small, fixed redshift and luminosity bins, respectively. 



Our results are fully consistent with the assumption 
that the X-ray variability mechanism and the accretion ef- 
ficiency are the same in all AGN, at all redshifts. In other 
words, our results strongly suggest that the X ray source 
operates in the same way in all AGN, at all redshifts. We 
find that, as the redshift increases, same luminosity AGN 
have smaller BH mass and higher accretion rates. At the 
same time, as the redshift increases, the rest frame light 
curve length decreases accordingly. For this reason the vari- 
ability amplitude (for the same luminosity AGN) increases 
slightly up to z~l, and then remains roughly constant with 
redshift. 

On the other hand, in the overall "variability 
amplitude redshift/luminosity" plots (shown in Fig. 2) the 
high luminosity objects are mainly those with the highest 
redshift as well. According to our results, their accretion 
rate should be the highest among the sources in our sam- 
ple. Consequently, they should also show a largo variability 
amplitude. However, they also have large BH masses, and 
their rest frame light curve length is also small. These two 
effects reduce significantly the observed variability ampli- 
tude, and can explain the global anti-correlation we observe 
between variability and redshift/luminosity. 

7.1. BH mass and accretion rate estimation 

To the extent that all AGN vary like the nearby ones (as our 
results suggest) , wc can provide BH mass and accretion rate 
estimates for the objects in our sample. Obviously, given the 
large scatter of the cr^xs values, such a measurement will 
not be very accurate if it were to be quoted individually for 
each object in our sample. The variability method can only 
provide average estimates of Mbh and riiEdd- 

Even in this case, these estimates can only be considered 
as indicative at this point. For example, the largest num- 
ber of sources (21) are in the "Low-z" bin. The riiEdd — 0.4 
model curve in the middle panel of Fig. 2 appears to agree 
well with the average [log((T^^g), log(L-x)] points in this bin, 
but this is hard to quantify. With just 3 points to com- 
pare with, we cannot perform a proper fit to judge the 
goodness of the model fit and provide confidence ranges for 
the best-fitting model parameter values. We need signifi- 
cantly larger samples of high-z AGN, in order to improve 
the accuracy of our results. 

Having these comments in mind, the model curves in 
Fig. 3 suggest that, on average, the Low (z^O.7), High 
(z-^l.S) and vHigh-z (z~2.9) Lockman Hole AGN accrete 
at ~ 25,40 and 50% of the Eddington limit. Our results 
are in agreement with those of McLure & Dunlop (2004) 
who also observe a slow evolution of the Eddington ratio 
from riiEdd ~ 0.15 at z~ 0.2 to riiEdd ~ 0.5 at z~ 2, for 
a large sample of quasars drawn from the SDSS catalogue. 
Although a few objects in Fig. 3 show variability amplitudes 
well above the average (by a factor up to 10 in some cases), 
these outliers are probably of stochastic nature (notice for 
example that the scatter of points at a given luminosity in 
the middle panel of Fig. 2 is of a similar magnitude). We 
conclude that, most probably, the majority of the sources 
in our sample are accreting at rates below the Eddington 
limit. This limit seems to be a relevant physical boundary 
to the AGN accretion rate out to redshift ~ 3. 

Regarding black hole masses, the model curves in Fig. 3 
suggest a BH mass range of 5 x 10^ — 2 x 10* Mq in the 
case of the Low-z objects, 1.5 x 10'' - 6.6 x 10* Mq for the 
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High-z objects, and 5 x lO'^ - 1.3 x lO'^ Mq for the AGN 
in the vHigh-z bin. Interestingly, even for these high-z and 
very luminous objects, the highest BH mass estimate is less 
than 3 x 10^ M©, i.e. the most massive BH measured dy- 
namically in the local Universe and the expected BH mass 
limit based on the known properties of early-type galax- 
ies and the locally observed correlation between bulge and 
black hole mass (McLure & Dunlop, 2004). 
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Appendix 

When dealing with variability, one of the main concerns 
is the possible contribution to the observed variations of 
systematic effects due to instrumental response and back- 
ground. Subtle variations of these instrumental quantities 
may remain uncorrected and can be difficult to spot. In or- 
der to verify that no systematic instrumental variations are 
introducing artifacts in the light curves we use in this work, 
we divided the light curve of each one of the V-objects (in 
the observer's frame) over its mean. Fig. 5 shows these "nor- 
malized" light curves. Most of the V-objects were mainly 
observed during the first five and the last three of the 10 
XMM-Newton observations listed in Table 2. This is due to 
the facts that a) we chose objects with Njata > 4 and b) 
the shift between the pointing coordinates of the October 
23 and 25 observations and the coordinates of the other 
eight observations is quite large. 

We then estimated the mean of all the normalized count 
rates within each observation (the "mean normalized count 
rate" of the XMM-Newton observations are listed in the 
second column of Table 3). If there are any significant bi- 
ases introduced by instrumental effects in any of the XMM- 
Newton observations, we would expect the respective mean 
normalized count rate to be significantly different than 
unity. It is only during the October 23observation that we 
observe the mean to be higher than unity, but even in this 
case this difference is significant at the 2.5(7 level only. This 
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Fig. 5. Top panel: The normalized light curves of all 
the "V" objects in our sample (open circles). Time is in 
Modified Julian Days. For clarity reasons, we have changed 
the x-coordinates randomly, by a small amount. Filled 
squares show the mean of all the points within each ob- 
servation (errors are also plotted, but are smaller than the 
symbol's size). Bottom panel: The average CCF between 
source V148 and the light curves of all the other "V" sources 
(in the observer's frame). 



result suggests that instrumental effects do not affect sig- 
nificantly the observed light curves. 

This conclusion is further strengthened when we con- 
sider the cross correlation between the light curves of the 
"V" sources in our sample. In the bottom panel of Fig. 5 we 
show the average cross-correlation function (CCF) between 
the V148 source and all the other V-sources (V148 is the 
brightest among the variable sources with Ndata = 8). To 
construct this average CCF we computed first the individ- 
ual CCFs between the V148 and the other sources, using 
light curves in the observer's frame, the DCF method of 
Edelson & Krolik (1988) and a lags bin of 2 days. We then 
computed the mean CCF at each lag where there were more 
than 10 individual CCF estimates. There is no obvious sig- 
nificant structure in the average CCF, which is rather flat, 
and consistent with zero at all lags. We reached the same 
result when we constructed the average CCFs using a few 
other "reference sources". We therefore conclude that any 
unaccounted instrumental effects should not contribute sig- 
nificantly to the variations we observe in the V-sources light 
curves. 

The third column of Table 3 lists the "root mean square 
(rms) variability amplitude", Fvar, of all the normalized 
count rates within each observation and its error (which ac- 
counts only for the effects of the flux measurement errors). 
Both the rms variability amplitude and its error have been 
estimated as in Vaughan et al. (2003). In effect, Fyar mea- 
sures the scatter of all the normalized count rates around 
their mean within each observation. In theory, we should 
expect that Fvar should remain roughly constant: we can 
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Table 3. The mean normalized count rate of all the "V" 
objects in each XMM-Newton observation, and the corre- 
sponding fractional root mean square variability amplitude, 



Obs Mean Fvar 

2002-10-15 1.10(±0.16) 0.51(±0.05) 

2002-10-17 0.98(±0.09) 0.44(±0.05) 

2002-10-19 1.11(±0.04) 0.15(±0.04) 

2002-10-21 0.87(±0.05) 0.25(±0.03) 

2002-10-23 1.23(±0.09) 0.35(±0.03) 

2002-11-27 0.99(±0.04) 0.25(±0.02) 

2002-12-04 0.94(±0.04) 0.24(±0.02) 

2002-12-06 0.99(±0.04) 0.24(±0.02) 



consider the flux measurements of the various V-sources in 
each frame as an ensemble of different measurements of the 
same object. If the variability process is stationary, then we 
would not expect the rms variability amplitude to change 
significantly. However, this is meant to hold "on average" , 
as Fvar measurements from individual light curves of an ob- 
ject can easily vary by a factor up to ^ 3 (sec for example 
the bottom two panels in Fig. 6 of Vauglian ct al., 2003). 

The values listed in Table 3 show that Fvar varies 
by a factor of ~ 2 — 3 from one observation to the 
other. It is difficult to understand how any instrumen- 
tal effects could introduce such a scatter in the F-^ar val- 
ues. However, even if this scatter in the Fvar values is 
the result of the red-noise character of the X-ray vari- 
ability in AGN, we considered the possibility it may af- 
fect the excess variance we measure for the faintest ob- 
jects in our sample. For that reason, we estimated the 
average Fvar of the 15 brightest and faintest objects in 
the V-samplc. More than 10 of them were detected dur- 
ing the October 17 and 23 observations (which show large 

^ var, total values). We find that Fvar, bright, Oct. 17 = 0.41 ±0.05 
and Fvar,bright,Oct.23 = 0.38 ± 0.05, while Fvar,faint,Oct.l7 = 

0.35±0.11 and Fvar, faint, Oct. 23 = 0.28±0.10. Given the fact 
that the fainter objects show smaller scatter around their 
mean, we are confident that instrumental effects do not in- 
troduce significant "artificial" variations in the light curves 
we use in this work. 



